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^ ■ Abstract 

oo 

The unit of quantum information is the qubit, a vector in a two-dimensional 
Hilbert space. On the other hand, quantum hardware often operates in two- 
K* ' dimensional subspaces of vector spaces of higher dimensionality. The presence 

Qs ' of higher quantum states may affect the accuracy of quantum information 

CS ■ processing. In this Letter we show how to cope with quantum leakage in 

devices based on small Josephson junctions. While the presence of higher 



o 



0\ \ charge states of the junction reduces the fidelity during gate operations we 



o\ 



demonstrate that errors can be minimized by appropriately designing and 



C^ ' operating the gates. 
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The most widely accepted paradigm of quantum computation describes quantum infor- 
mation processing in terms of quantum gates whose input and output are two-state quantum 
systems called qubits ^. Quantum Computation (QC) is performed by means of a control- 
lable unitary evolution of the qubits 0. Due to the intrinsic quantum parallelism, problems 
which are intractable on classical computers can be solved efficiently by using quantum 
algorithms. Probably the most striking example is the factorization of large numbers [||. 

Parallel to the developement of the theory of quantum information there has been an 
increasing interest in finding physical systems where quantum computation could be imple- 
mented. In an (almost) ideal situation one should identify a suitable set of two-level systems 
(sufficiently decoupled from any source of decoherence [Q) with some controllable couplings 
among them needed to realize single qubit and two-qubit operations. These requirements 
are sufficient to implement any computational task 0. Various physical systems have been 
suggested for the implementation of quantum algorithms, e.g. ions traps 0, QED cavi- 
ties [0 and NMR [^. The quest for large scale integrability and flexibility in the design 
has very recently stimulated an increasing interest in the field of nanostructures. Up to 
now promising proposals are based on small-capacitance Josephson junctions 0-0, cou- 
pled quantum dots |[T^JT5[| and phosphorus dopants in silicon crystals [0 . The experiments 
on the superposition of charge states in Josephson junctions |T^JT^ and the recent achieve- 



ments in controlling the coherent evolution of quantum states in a Cooper pair box |T^ 
render superconducting nanocircuits interesting candidates to implement solid state quan- 
tum computers. 

Physical realizations of QC are never completely decoupled from the environment. Since 
decoherence will ultimately limit the performance of a quantum computer a lot of attention 
is being devoted to this problem. Besides decoherence, for each proposed scheme a detailed 
analysis of the errors induced by the gate operations themselves is crucial in order to assess 



their reliability and the feasibility of fault-tolerant quantum computation ||2^J21|. Errors 
may occur due to a variety of reasons. An obvious example are fluctuations in the control 
parameters of the gate which act as a random noise and thus affects the unitarity of the 
time evolution. Alternatively gate operations can change the coupling of the qubits to the 
environment (even if this coupling is negligible during storage periods) thereby enhancing 
decoherence. All these error sources can be analyzed by properly modelling the qubit- 
environment coupling. However, there are errors which are not due to (or cannot be described 
in terms of) the action of an external environment. Much rather they are inherent in the 
design of the gate. 

In this Letter we consider one (intrinsic) source of error in gate operations which is 
common to several of the proposed solid state implementations, the quantum leakage. It 
occurs when the computational space is a subspace of a larger Hilbert space. This is the 
case e.g. when the information is encoded in trapped ions or in charge (or flux) states 
of devices based on Josephson junctions (or SQUIDS). We start by introducing a general 
scheme to characterize the leakage and then we focus on devices based on small-capacitance 
tunnel junctions. 

Our analysis applies to the situation illustrated in Fig. |l|. The two low-energy states 
constitute the computational Hilbert space. The system, however, can leak out to the 
higher states. If the energy difference between the low-lying and the excited states is large 
compared to the other energy scales of the system (as in Refs. [0-0) the probability to leak 



out is small. One might wonder whether it is necessary to discuss this effect at all. As we 
will see the consequences of leakage are more severe than a simple estimate of energy scales 
might suggest. The presence of states outside the computational space modifies the time 
evolution of the qubit states compared to the idealized design. 

The ideal unitary gate operation Uj is obtained by switching on a suitable Hamiltonian 
Hj which couples the desired computational states in a controlled way for a time to- By 
choosing to one can implement the desidered gate operation. In reality, however, the dy- 
namics of the system is governed by a unitary operator Ur which acts on the full Hilbert 
space. Since information is being processed within the computational subspace the output 
is related to the input state via the map nf/R(t)n, where 11 is the projection operator on 
the computational space. One is interested in optimizing the real gate operation in order to 
get as close as possible to the ideal [//. In general the "best" operation may require a time 
t 7^ to as all the system eigenenergies are modified by the states outside the computational 
subspace. Therefore we use the time t as parameter to optimize the given computational 
step. We characterize the performance of real gates by the fidelity JF and the probability of 
leakage £(t) defined as 

T=i- ^min^t}llf^/(^o) - nf/K(t)n|| (1) 



Cit) = 1 - mm4i;\Ulit)nURitm (2) 

In Eq. ([^) we make use of the operator norm defined as \\D\\ = Sup^\D\ 



Sup^J {tplD'^Dlip) over the vectors {{ip) : {ip\tp) = 1} of the computational subspace. This 
definition implies that \\D\\ = \/\m where Am is the biggest eigenvalue of D'^D. As in 
the case of the minimal fidelity [^] this definition gives estimates for the worst case. The 
definition given in Eq. (|1|) can therefore be regarded as a prescription how to optimize the 
gate design (note that the fidelity defined in Eq.(|I|) does not depend on the time t). 

As mentioned before the existence of states other than the computational ones has two 
main consequences on the qubit dynamics. There is a nonzero probability of leakage, mea- 
sured by £(t), and a modification of the eigenenergies and eigenstates of the real system. 
The latter effect turns out to be an important source of gate errors. 

In order to study the phenomena related to leakage quantitatively we apply Eqs. (P, 
(0) to Josephson junction qubits in the charge regime as proposed in Refs. [p|,p^. A similar 
analysis can be carried out, with appropriate changes of paramaters, for all other cases where 
leakage is present. In Refs. [§,|r^ the qubit is implemented using nanocircuits of Josephson 
junctions. The corresponding Hamiltonian for one and two qubit operations can be written 
as 

Hr= Yl [Echini - rix^if - Ej cos (pi + El (sin 0i + sin ^2)^ (3) 

j=l,2 

In the first term E^^ is the charging energy. The second and the third term represent the 
Josephson tunneling (associated with the energy Ej) and the inductive coupling of strength 
El [^ which bring about single and two qubit operations possible. Both Ej and E^ 



are assumed to be much smaller than the charging energy. The offset charge n^i can be 



controlled by an external gate voltage. The phases 0j and the number of Cooper pairs rzj 
are canonically conjugate variables [02,%] = i Sij [24 . 



At temperatures much lower than the charging energy, for Ux^i ~ 1/2 the two charge 
states rii = 0,1 are nearly degenerate. They represent the states |0), |1) of the qubit (see 
Fig. 0). In the computational Hilbert space the ideal evolution of the system is governed by 
the Hamiltonian 
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Hi = ^ AEch,i(Jz,i 7r(^x,i 7Z-C^y,lC^y,2 (4) 
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where a are Pauli matrices and AEch,i = Ech{nx,i — 1/2)- The different time evolution 
due to Hn and to Hj causes an error in the gate operation. We note that leakage is also 
present during idle periods of the gates. However, here we only discuss the errors during 
gate operations. 

One-bit gate {E^ = 0). Single-qubit gate operations can be implemented, e.g., by sud- 
denly switching the offset charge to the degeneracy point n^; = 1/2 where the charge states 
|0) and |1) are strongly mixed by the Josephson coupling ||2^. Whereas in the ideal setup 
this coupling mixes only the states |0) and |1), in the real qubit all charge states are involved. 

The evolution in the computational subspace for a time interval t is described by the 
operator {h = 1) 

nun{t)u = ^e-*^"*n|<i'„)(<i'„|n (5) 

where 11 = |0)(0| + |1)(1| is the projector on the computational subspace and |$n) are the 
eigenstates with energies En of the Hamiltonian H^ (here |$„) can be expressed in terms of 
Mathieu functions). 

By evaluating the leakage according to Eq. (^ we obtain 

/:(t) = 1 - min± I J2 (±r(0|$„)($Jm)e-^^"*|' 

n;m,=0,l 
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[1 - min± cos(2Ech ± Ej/2)t] . (6) 

The order of magnitude {Ej/ Echf can be understood immediately by regarding the coupling 
to higher charge states as a perturbation to the ideal system of Eq. @). 

The fidelity has to be hmited by the leakage since it describes the length of the projection 
of the true state at time t onto the ideal state at to- There is another effect contributing to 
the loss of fidelity: the presence of higher charge states renormalizes the energy eigenvalues 
thus leading to a frequency mismatch between ideal and real time evolution. However, due 
to the symmetry of the system and the fact that Ej is the only coupling energy to the 
states outside the computational subspace there is a simple way to cure this problem. Let 
us consider a 7r-rotation. The optimal gate is obtained by changing the operation time to 
tp = -K I AE where IS.E is the energy splitting between the two lowest eigenstates (as opposed 
to the time to = '^ jEj in the ideal system). The value of the fidelity is then given by 
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l--^^2 + 2\sm{2nE,^/Ej)\ (7) 



We mention that the error accumulates hnearly with the number of operations. For typical 
parameters of Josephson junctions Ej/Ech ~ 0.02 one finds that after about 10^ operations 
the loss of fidelity becomes of order unity. 

Two-bit gate {E^ ^ 0). Among the many possibilities for the elementary two-qubit op- 
eration, choosing a particular one may be a non-trivial step in the course of implementing 
quantum hardware. Due to universality of quantum computation [^] one is free to use any 
generic 4x4 unitary matrix as a two-qubit gate. From our point of view a choice is op- 
timal if it avoids errors stemming from a discrepancy of the ideal gate and the way of its 
implementation. Therefore, in the following we assume that the Hamiltonian as introduced 
in Eq. (§) 
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-Ej/2 -El/2 -Ej/2 
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(8) 



generates the ideal two-bit gate. In Eq. (^ we have used the basis {|00), |01), |10), |11)} 
(which is obtained as the direct product of the states introduced previously). The typical 
scale for the operation time to is on the order of I/El P,P • 




In complete analogy with the one-bit gate we find that the leakage is of the same order 
for the two-qubit operation: 

C oc max 

(the numerical coefficient is larger than in the one-bit case because there are more charge 
states outside the computational subspace directly coupled to the qubit states either by El 
or Ej). 

The situation for the fidelity, however, is different. In order to estimate J-" we consider a 
perturbative expansion of D'^D where D = t//(to) — IlULi{t)Il up to second order in Ej/E^h, 
EL/Ech- The eigenvalues of this matrix have the form 2 — 2cos(-E'„t — E^Hq) + 2nd order 
terms (here En and E^^ are the eigenvalues of Hr and Hj, respectively). It turns out that 
due to the presence of several energy scales the frequency mismatch between real and ideal 
time evolution cannot be compensated for by adjusting the operation period. The leading 
terms of the fidelity can be written as 

^^l-l/'g ^J +b^] , (9) 

2 \ ElEcIi EchJ 

where a and b are coefficients which depend on the particular choice of rix^i and to- In Fig- H 
we show the numerical results for rix = 1/4 and to = t^/ El- The loss of fidelity (the term 
in paranthesis in Eq. (P)) is proportional to to. The maximum (the best operation one can 
achieve) scales linearly with Ej/Ech- This should be contrasted with the one-bit case where 
it scales quadratically. 

We mention that we have chosen the definitions for the leakage and the fidelity describing 
the "worst case" in order to avoid a dependence of the discussion on the preparation of the 
initial state. One could wonder whether the "generic case" is much more robust with respect 



to leakage. It is easy to convince oneself by checking various choices of initial states that 
the loss of fidelity is indeed on the order of the worst case estimates. 

In conclusion, starting from given gate operations we have discussed their optimal im- 
plementation in real systems. We have shown that leakage limits the number of operations 
which can be performed reliably both for one and two qubit gates. For one-bit gates one 
can correct leakage errors by changing the operation time. We have pointed out that with 
respect to fidelity it may be appropriate to choose the elementary two-qubit gate as it is 
determined by the implementation. Fig. ^ shows the central result of this work: although 
leakage causes an inevitable loss of fidelity for two-qubit operations, this loss can be mini- 
mized by an appropriate choice of the device parameters. 

Finally we mention that one can speculate about correction procedures for errors caused 
by leakage. It should be possible to check during the computation whether leakage has 
occured. This should be done by measuring the system only if it is outside the computational 
subspace. One can imagine to realize a low sensitivity SET transistor which is able to 
measure the system only if the charge is outside a specified window. 
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The coupling E^ is also proportional on the Josephson coupling energys. For our pour- 
poses we can assume that is an independent energy scale in the problem. 
The model ignores quasiparticle tunneling which is suppressed if the temperature is 
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do not examine the various possibilities to realize one and two-qbit gates. The qualitative 
results do not change although it may be more convenient, as far as leakage is concerned, 
to choose a particular scheme. 
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FIG. 1. Schematic view of a qubit with leakage. The two low energy states constitute the 
computational Hilbert space. The system however evolves with the opertor Ur and therefore can 
leak out to the higher excited states. In the case of Josephson junction leakage is due to the 
Josephson tunneling to high Cooper pair charge states. In the case of two qubit operations the 
computational space is spanned by the states {|00), |01), |10), |11)} and the coupling with the higher 
charge states is due both to Ej and to El 
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FIG. 2. Fidelity for a two-qubit gate as a function of El for different values of Ej. 



